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Abstract 

We investigate a perturbation of a scalar field model (called here the 
signum-Gordon model) with the potential V(f) = \f\. The perturbation 
generalizes the signum-Gordon model to the signum-Klein-Gordon model 
i.e. to the case V(f) = \f\ — \\f 2 , where A is a small parameter. Such 
a generalization breaks the scaling symmetry of the signum-Gordon model. 
In this paper we concentrate on solutions for self-similar initial data. Such 
data are particulary useful for identification of the effects caused by the term 
that breaks the scaling symmetry. We have found that the behaviour of the 
solutions is quite interesting - they escape and return periodically to the self- 
similar initial data. 
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1 Introduction 



The present paper refers to recently investigated scalar field models with V-shaped 
potential dTJ-OJ . Such potentials have a common feature - left- and right hand 
derivatives are different from zero at the minimum. Mentioned models have a 
well-justified physical origin despite the fact that they seem to be a little bit ex- 
otic from a mathematical viewpoint. Moreover, some physical systems described 
by scalar field models with V-shaped potentials are easy to built (e.g. chain of 
pendulums impacting on a rectilinear bar). Unfortunately, such models have a 
very unpleasant mathematical feature - a typical solution consists of many (some- 
times infinitely many) partial solutions. The partial solutions are matched at some 
points. The matching procedure is mostly onerous. This is probably the reason 
why literature is poor in results for field theoretic models with V-shaped potential. 

It turns out that the behaviour of the field close to minimum strongly depends 
on a 'shape' of the potential [4J. In particular, for V-shaped potential a field ap- 
proaches exactly its vacuum value at finite distance (a parabolic approach). This 
fact has a profound significance - kinks have no exponential tails! Such kinks are 
called compactons because their supports are compact. The compactons consid- 
ered in our models are topological, so they are qualitatively different from e.g. 
well-known compactons in the modified KdV model Q-flTl- Recently, the topo- 
logical compactons have been also obtained in models with nonstandard kinetic 
terms as so-called fc-defects [8J. It is important to notice that there are other (non- 
topological) compact solutions in the s-G model [9|. Because of their properties 
they are called oscillons. 

The second characteristic property of the mentioned models is a scaling sym- 
metry, see O. This symmetry means that if a function f(x, t) is a solution of a 
field equation, then new function defined as f u (x, t) = u 2 f(x/u, tjv) is a solution 
as well, v is here a positive constant. A presence of the scaling symmetry in the 
model suggests existence of solutions that are invariant with respect to the scaling 
transformation (so-called self-similar solutions). Such solutions have been ob- 
tained in the s-G model. A complete list of solutions for the self-similar initial data 
is presented in ifTOl . For models with the potential V(f) = afQ(-f) + bfQ(f), 
where a, b are constant parameters and is the well-known Heaviside step func- 
tion, the scaling symmetry is exact whereas for most models with V-shaped po- 
tential the symmetry is only approximated. Note that the signum-Gordon model 
(s-G) can be obtained as a particular case, i.e. by setting a = — 1 and b = 1. In a 
group of models with symmetric V-shaped potential, the s-G model is the simplest 
one. In this paper we study just symmetric potentials. 
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The aim of the present work is the analysis of the perturbed s-G model, where 
for simplicity reasons the specific perturbation is chosen in the simplest, nontrivial 
form. Namely, we add the quadratic term. In spite of its simplicity such a gen- 
eralization of the s-G model allows to face several important problems. The first 
one is breaking of the scaling symmetry. Among the physical systems there are 
fewer of them with an exact scaling symmetry. There are always fluctuations in a 
typical physical system that interacts with its environment. The fluctuations mod- 
ify an effective potential and break the exact scaling symmetry. In this physical 
context, it is clear that investigation of the perturbed field theoretic models with 
V-shaped potentials is an important issue. It allows for better understanding of 
dynamics of compactons in the systems with the broken scaling symmetry. In our 
paper we analyse the perturbed potential V(f) — \f\ — \\f 2 , where A is a small 
parameter i.e. |A| <C 1. The second important problem, which is in general very 
difficult for systems with non-differentiable potentials, is a stability analysis of 
solutions. Our investigations are some kind of structural stability analysis. Such 
analysis is important for compact kinks as well as compact oscillons. 

However, the models with V-shaped potentials are interesiting from the math- 
ematical point of view, they have also some properties that allow to think about 
possible applications to condense matter physics and cosmology as well. In the 
cosmological context, the most interesiting seems to be the fact that for models 
with potentials sharp at its minima the terms that come from a gradient of the 
potential dominate the field dynamics close to the minimum. For instance, in the 
s-G model the term ^ = sign/ remains finite arbitrary close to the minimum. 
This is in total oposition to the behaviour of e.g. 4 theory, where gradient of 
the potential vanishes close to the minimum. Because of this, small perturbations 
propagate easily within the topological compactons or other nontopological field 
configurations like, e.g. mentioned oscillons, whereas outside of them the prop- 
agation encounters on resistance. Moreover, an absence of linear perturbations 
around the V-shaped minimum is a basic feature of our models. It entails auto- 
matically that the linear perturbations can propagate only at a defect background. 
This effect, characteriscic for models with V-shaped potentials, is similar to be- 
haviour of A;-fields that play a prominent role in cosmology (see flTU). 

In our calculation we concentrate on differences between solutions in the s-G 
model and solutions in the signum-Klein-Gordon (s-K-G) model (the perturbed 
model). Applying the same initial data for solutions in both these models we can 
analyse the differences between their solutions as a pure effect caused by the term 
\ \f 2 . In the case when initial data are self-similar (parabolic) the solutions in the 
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s-G model have especially simple form. We apply the self-similar data just for 
this reason. 

Our paper is organized as follows. In the section 2 we introduce the signum- 
Klein-Gordon (s-K-G) model and give a general method of calculation of partial 
solutions that can be obtained directly from the self- similar initial data. Unfor- 
tunately, they are insufficient to construct a solution valid for each x and t > t , 
where to is an initial moment. Section 3 is devoted to a study of a solution for a 
specific self-similar initial data. Focusing on a specific initial data enables us to 
calculate all partial solutions that (when matched together) cover the whole range 
of variable x. In the last section we summarize our results and emphasize effects 
that stem from the term \\f 2 . 

2 Initial problem for the generalized model 

2.1 The signum-Klein-Gordon model and its partial solutions 

The s-K-G model for the scalar field f(x, t) in 1+1 dimensions has the Lagrangian 
L = \{d t ff~ l -{d x ff-V{f\ (1) 
where the potential V(f) is given by the formula 

V(f) = |/| - \\f 2 . (2) 

Euler-Lagrange equation that corresponds to Lagrangian (OQ) takes the following 
form 

(d 2 - dl)f + sign/ - A/ = 0. (3) 

The sign of parameter A has a crucial meaning for the behaviour of the field 
f(x, t). The potentials V(f) for negative and positive values of parameter A are 
qualitatively different Q. The case A = gives the s-G model which has been 
discussed in our previous papers, (see e.g. fll, GO and ifTOl ). In this paper we are 
interested in the case A < 0, because for A > the potential V(f) is not bounded 
from below. Nevertheless, a perturbative method presented in the following sub- 
section involves both cases of sign A. In order to distinguish between different 
kind of solutions, we use symbol / for solutions in the s-K-G model (A ^ 0) and 
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symbol (f) f° r solutions in the s-G model. The partial self-similar solutions are 
given by the formula 

w= tttjtfcsa, ,4) 

2 Vk^Vk - 1 

where x £ fufc_it, i>fct] and fc = 1, 2, The parameters Vk are velocities of zeros 

of polynomials. They are determined from matching conditions. For more details 
see [10]. The partial solutions obey the relation sign0 fc = (— l) fe+1 . By analogy, 
we define partial solutions in the model with A ^ 0. They obey the equation 

(fi? - %)f k (x,t) - (-l) fc - Xf k (x,t) = 0, (5) 

where sign/ fe = (-l) fc+1 . For |A| < 1 the potential V(f) = \ f\ - \\f 2 can be 
interpreted as a perturbed potential V(<f>) = \(f)\. In this case we say that the exact 
scaling symmetry is violated or that the generalized model has an approximate 
scaling symmetry when \ f\ ^> A/ 2 . 



2.2 Self-similar initial data and partial solutions 

This paper is devoted to investigation which are the effects caused by the term A/ 
in the s-K-G equation. It can be achieved by comparison two solutions for the 
same initial data: the first one that is a solution in the s-G model and the second 
one that comes from the s-K-G model. The differences between them are a direct 
consequence of the term that breaks the scaling symmetry. From practical reasons 
we investigate some characteristic points of solutions, i.e. trajectories of its zeros. 

It turns out that explicit formulae for the solutions are not always available - 
this problem strongly depends on initial data. It has been shown, see IfTOl . that 
solutions in the s-G model for the self-similar (parabolic) initial data are given 
by explicit formulae. For this reason, the self-similar initial data are more useful 
for our purposes than other, more general, initial data. In fact, any self-similar 
solution <fi at the moment t = t is suitable for our purposes and can be applied as 
an initial data. Therefore, we assume following initial data for the partial solutions 

fk(x,to) = (j) k (x,t ), dtfk(x,t)\ t =t = d t (j)k(x,t)\t=t Q - (6) 



2.3 The perturbative method 

The method presented in the current subsection allows us to obtain the partial 
solutions fk(x, t) directly from the initial data ©. We call them the partial so- 
lutions of the first kind. It turns out that such partial solutions are insufficient. A 



5 



complete solution f(x, t) is consisted of some other partial solutions as well. This 
inconvenience appears also for, e.g. the s-G equation in the case when initial data 
have the form of piecewise smooth functions matched up at some points. In most 
cases, such matching points are origins of new partial solutions. At the initial mo- 
ment t = t mentioned partial solutions are shrunk to single points but for t > t 
their supports expand (the partial solutions of the second kind). The s-K-G model 
has analogical partial solutions (the first and the second kind) - we discuss their 
properties in the further part of this paper. The partial solutions discussed in this 
paragraph are the first kind ones. They cover whole range of axis x at t = t (they 
obey © where 4> k have this property) and their supports shrink for t > to- 

Let us assume that the solutions f k (x,t), which depend on parameter A, are 
represented in the form of a power series 

oo 

f k (x,t) = J2fkn(x,t)\ n , (7) 
n=0 

where A< 1. After plugging series © into equation © we get a set of equations 

f (d 2 -d 2 x )f k0 (x,t) = (-l) k n = 0, 

\ (d?-d 2 x )f kn (x,t) = fkn.^t) 71=1,2,3,.... W 

Each of them takes the form of wave equation with a source. We can integrate 
them using new variables £ = |(x + t), rj — \{x — t). The result of integration in 
the original variables reads 

f k0 (x, t) = F (x + t) + G (x -t)- t^(x 2 _ t 2 ), (9) 



x-\-t x — t 

fkn(x,t) = F n (x + t) + G n (x-t) - / 2 da 2 d/3 jfw-i(ar + /3, a - P), 

Jo Jo 

(10) 

where F(x + 1) and G(x — t) are arbitrary functions. They can be calculated from 
the following initial conditions for the partial solutions: 

fko(x,t ) = <pk(x,t ), d t f k0 (x,t)\ t=t0 = d t (p k (x,t)\ t=to , (11) 
fkn(x,t ) = 0, d t f kn (x,t)\ t=t0 = 0. (12) 
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Conditions (fTTT) and (PT21) stem from the initial conditions ©. The solution © 
obeys the first of the equations ([8]) (i.e. the s-G equation) and the self-similar 
initial data (fTTT) . so it coincides with fa, 

f k o{x,t) = fa{x,t). (13) 

A direct calculation confirms this result. It means that the self-similar solutions 
(fT3l are zero-order approximation for the solutions ©. In order to get higher- 
order approximations we have to find the functions F n (x + t) and G n (x — t). 
Differentiating first of equations (PT21) with respect to x, combining with the second 
one and shifting arguments, we obtain equations 

D+f k n(x, t) \ x=s _ t0t t=tQ = 0, D_f kn (x, t) \ x=w+tQt t=to = 0, (14) 

where D± = \{d x ± d t ). Equations (fl4l) can be rewritten in the form 



x-\-t x — t 

F' n {s) = D + [ 2 da [ 2 d{3 hn-M +P,a-P) 
Jo Jo 



(15) 



X=S — tQ, t=to 



x-\-t x — t 

G' n {w) = D_ f~ da f~ d{3 hn.^a +/3,a-/3) 
Jo Jo 



(16) 

X=W+to, t=tg 



where formula (TTOl) has been applied. The sum of constants that comes from 
integration of the expressions (fT5l) and (fT6l) is fixed by the first of conditions (TT2l) . 

In the first step we calculate the function f k i(x,t) from formula (flOl) . Then 
we can continue the procedure in order to obtain /^(z, t). In principle, this pro- 
cedure can be repeated infinitely many times giving expressions for all functions 
fkn(x, t). In fact, we are able to obtain only few functions f k n(x, t) because the 
calculations became quickly too complicated. Fortunately, functions f kn (x, t) ob- 
tained for several, the lowest values of n enable us to guess a general formula for 
arbitrary n. This formula has the form 

f kn (x, t) = (-i) k i~ t ° ) 2 An: * +Bnx ~t Cn , 



where 



(2n + 2)\ ffc_iffc 

A n = (2n+l)(n + l), 
B n = ~{v k -i + v k )(n + l)(t + 2nt ), 
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C n = (n + v k _iv k )t 2 + n(l + (2n + l)^^)^- 

One can check that formula (fTTT ), which was originally found for n = 1, 2, . . ., is 
also true for n = - in this case it gives (IT3T ). 

It turns out, and this is a big surprise, that the series © can be summed up 
giving as a result 



fk(x,t) 



+ 



2(v k - 1 v k 



Mcosh (p(t-io)) + N 



sinh (p(t- t )) 



+ 



[cosh (/>(*- to)) - 1], 



(18) 



for A = p 2 > 0, and 

/*(<M) = 



2(f fc „iw fc - 1) 
("l) fc 



Mcos (cr(t - t )) + iV 



sin (a(t - 1 )) 



(7 



£7 



1 - cos(cr(t -t ))], 



(19) 



for A 



-a 2 < 0. The coefficients M and iV read 

M = (x - v k -ito)(x - v k t ), 

N = t-t + 2v k - 1 v k t - (ffe-i + v k )x. 

The partial solutions f k satisfy the relation sign/^ = (— l) fc+1 . 

It is worth emphasizing that the partial solutions f k in the model with the 
explicitly broken scaling symmetry, still have a quadratic dependence on variable 
x. What changes is the time dependence. This is the first important result obtained 
with the help of the perturbative method. 

Knowing this result we can, of course, propose a posteriori a proper Ansatz 



f k (x,t) = a{t)x 2 + b(t)x + c(t). 



(20) 



Plugging the Ansatz (|20l ) into eq. © we get the set of ordinary differential equa- 
tions for the coefficients a(t), b(t) and c(t) 



d 2 a 



Xa = 0, 



Xb = 0, 



Xc = 2a+{-l) k . 



The constants that come from integration of these equations are fixed by the con- 
dition ©. Solving these equations we recover formulae (TTSl and <fT9l depending 
on sign of the parameter A. 
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3 The solution for a specific initial data 



In this section we present a solution in the s-K-G model for a specific self-similar 
initial data. We concentrate on the case when the perturbative parameter A is 
small |A| <C 1 and negative A = —a 2 . The main effort is focused on the partial 
solutions of the second kind. They can be obtained as the solutions of boundary 
problem. The boundary conditions are given by the partial solutions of the first 
kind at points of contact of their supports and supports of the partial solutions 
of the second kind. The most serious obstacle in the s-K-G model is that we do 
not know a general formula for the partial solutions. For instance, in the s-G 
model such formula consists of two arbitrary functions and terms ±t 2 /2, ±x 2 /2 
or their combinations. For this reason, we restrict our study to a specific self- 
similar initial data and construct an approximated solution. Nevertheless, it turns 
out that investigation of such a specified case gives valuable information as well. 
One of the most interesting results presented in the current section is a discovery 
of periodicity in time for such solution in the s-K-G model. 

3.1 The positive partial solution /+(#, t) 

Among the self- similar initial data, the simplest one reads 



where Q(x) is the Heaviside step function. The self-similar solution cf>(x, t) in the 
s-G model for data (|2~TT) is consisted of two partial solutions matched up at the 
light cone x = t, see also [10J. The solution (f>(x, t) has a very simple form 



This formula corresponds to © for k — 1, v — — 1 and Vi — 1. The snapshot of 
the solution <f){x, t) is presented in Figffl 

In the further part of this paper we concentrate on the solution f(x, t) in the 
s-K-G model. The positive partial solution f+(x, t) obeys the equation 



(f)(x,0) = -x 2 B(x), d t (j){x,t)\ t= o = 0, 



(21) 




(22) 



(d 2 -dl)U + a 2 f + + l = 0, 



(23) 



and the initial conditions 



f + (x,0) = ^x 2 , 



d t f+(x,t)\t = o = 0. 
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x = t 



Figure 1 : The self-similar solution for initial data (121) . The dashed line represents 
an initial configuration of the field <p. 



It takes the form 

f+(x, t) = - cos(at)x 2 + 4- sin(at) - ^(1 - cos(trt)). (24) 
4 4<T cH 

The partial solution (1241) can be obtained directly from <fl~9l ) for = 1 and t = 0. 
The formula (1241) holds for x > The trajectory of zero x\(t), 

^i(^) = "A/ - ^tan( C rt) - 1 (25) 

a y cos(crt) 4 

is the solution of the equation f + (xi,t) = 0. Except for the point t = 0, the 
function xi(t) obeys inequality xi(t) > t, what means that zero of f+(x, t) moves 
faster than its counterpart (x(t) = t) in the s-G model. Moreover, we can see that 
the velocity of zero xi(t) depends on variable t. Let us remind that velocities Vk of 
zeros of the self-similar solutions are constant, see formula ©. It means that the 
term \f in the eq. © is responsible for non monotonous expansion or contraction 
of the supports of the partial solutions. 

A series expansion of the formula (|25l) for small t 

x 1 (t)=t + W + ^tW + 0(e) (26) 

gives valuable information as well. We can see from (|26l ) that the zero Xi(t) 
moves with the acceleration x\ (t) which is proportional to a 2 provided that i < 1. 
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Figure 2: The trajectory xi(t) for A = —0.1. 



This observation has a practical meaning because it enables us to calculate the 
parameter A from experimental data. The curve X\(t) is presented in the Fig. 
[21 We can see a very good agreement between the analytical calculation and the 
numerical data. The function (12~5l) goes to infinity for t — > t* = ^- . It means that 
the solution f+(x, t) is valid for t < t*. In our numerical calculation (a 2 = 0.1) 
the characteristic time t* ~ 4.9673. The leading behaviour of xi(t) close to t* is 
given by the first term of the expansion 

1 / 16 — 2tt 3 / t* -t 1/2 . 

3.2 The negative partial solution f_(x,t) - some general re- 
marks 

In this and two further paragraphes we present the partial solution of the second 
kind f-(x, t). Our solution is only approximated and holds for times not longer 
than t ~ 2. Such a partial solution can be obtained as the solution of boundary 
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problem because at the initial time t = its support is a single point and it is 
located at x = 0. For later times t > the support expands to infinite size. This 
behaviour has been observed in our numerical calculation and it is suggested by 
the fact that xi(t) tends to infinity for t — ► The solution f-(x, t) has to obey 
the following boundary conditions 

/_(ari,t) = 0, d x f-(x,t)\ x=Xl = d x f + (x,t)\ x=Xl , (27) 

f-{x ,t) = 0, d x f4x,t)\ x=X0 = 0. (28) 

Conditions (1271) and (|28l) are derived from the field equation © and they mean 
that the partial solutions are matched so that the solution f(x, t) is smooth at xi(t) 
and xo(t). At xo(t) the partial solution f-(x, t) is matched with the trivial partial 
solution fo(x, t) = 0. Whereas the point x±(t) is given by the formula (|25l) . the 
second zero of f(x, t), i.e xo(t), is not known yet. In Section 3.4 we show how to 
obtain an approximated formula for x (t). 





f(x,t) 




/ 4 




u 










x = t 


x l / 




X 









==— 1 


/ fS 







Figure 3: The solution f(x, t) for initial data (I2T1) and times t < j- The dashed 
line represents an initial configuration of the field /, 

It turns out that the partial solution that obeys (|27l) does not obey (|28l) and vice 
versa. In order to get rid of this inconvenience we split the solution f-(x, t) into 
two pieces f^(x, t) and f R (x, t). They are matched at the light cone x = t. Such 
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split is sufficient to obtain /_ (x, t) that obeys (1271 ) and (1281) simultaneously. The 
solution f(x, t) consists of the following partial solutions 



f(x,t) 



for 

f^(x,t) for 

f*(x,t) for 

f+(x,t) for 



X < x (t), 

Xo(t) < X < t, 
t < X < Xi(t), 
X > Xi(t), 



where t < The snapshot of f(x, t) is depicted in Fig. [3]- compare it to the 
solution in Fig. \T\ 



3.3 The partial solution f R (x, t) 

It has been already mentioned at the begining of Section 3 that a general formula 
for the partial solutions in the s-K-G model is not known. This is the most serious 
obstacle in our investigations. Therefore, we search for the approximated partial 
solution f-(x, t). The approximation of eq. © is obtained by replacing the term 

cr 2 /(x, t) by the term a 2 (f)(x, t), what gives us 

(d 2 - d 2 x )Ux,t) + a 2 0(x,t) -1 = 0. (29) 

Such a modification is valid only for small times. The partial solution f^(x, t) of 
eq. (|29l (an approximate solution of eq. © for A = — a 2 ) at t < x < x 1 (t) has 
the form 

/*(*, t) = F R {x +t) + G R (x -t)- -(x 2 - t 2 ) + ^-{x 2 - t 2 ) 2 . (30) 

4 64 

In accordance with (1271) . this solution is matched to f+(x, t) at x — x\(t). It turns 
out that exact formulas for Fn(x+t) and Gr(x— t) can not be achieved because we 
need the inverse functions of x\ (t) ±t, where x\ (t) is given by (|25l) . Nevertheless, 
we can expand the expressions x\ (t) ±t in power series and then invert these series 
up to an arbitrary term. This is why we concentrate on series expansions of the 
partial solutions. In our further calculations we use the perturbation parameter a 
as an expansion parameter. The partial solutions represented by finite series (i.e. 
approximated partial solutions) obey the matching conditions (l27l) and (1281) up to 
some range of a. In order to find this range we start from series expansion of the 
formula (1251) for a 1, i.e., 

x 1 (t)=t+itV + ^V + 0(a 6 ). (31) 
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Note that the expansion (|3TT ) has the same form as the expansion (|26l ) for small 
times t. A leading term of the expression x\ (t) — t is proportional to a 2 . The pow- 
ers of the expression x\ (t) — t appear in ((271) because the partial solution f^(x, t), 
which is given by formula (1301) . includes terms proportional to x— t and a 2 (x— t) 2 . 
In order to take into consideration contributions from all terms in (1301) . especially 
from o 2 {x — t) 2 , we need terms proportional to a 6 at least. This is, naively, an 
accuracy of f^(x, t). The real accuracy is lower. The direct calculations allow us 
to obtain f^(x, t) only up to terms proportional to <x 4 because solutions of (1271) . 
i.e. F' R and G' R have such accuracy, see formulae below. The prime ' stands for 
derivatives with respect to whole arguments of Fr and Gr. 

In the first step we expand formulae (1241) . (1301 ) and their derivatives with re- 
spect to variable x at x = x\(t), where x\{t) is given by the formula (13TI ). The 
result has the following form 



/+(xi(t),t) =0(a 6 ), 



1 37 

/*(x lf t) = ^(xx + 1) + G R {x x -t)- -tV - — tV + 0(<x 6 



11 13 

<9*/*(x, t)!™ = F' R (xi + G'Rix, ~t)--t- -tV - — tV + 0(a 6 ). 

Plugging three last formulae into conditions (1271) we obtain two equations that 
contain F R , Gr, F' r and G' R . Then, we differentiate the first of these equations, 
i.e. f_(xx(t),t) = with respect to variable t. In the next step we solve obtained 
equations with respect to F' R and G' R , which gives 

F'r(xi +t) = - Va 2 - ^LtV + 0(<x 6 ), (32) 



1 67 

G' R ( Xl -t)=t+ -/a 2 + — tV + 0(a e ). (33) 
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Eq. (1321) can be integrated with the help of new variable s = xi(t) + t, where 
Xi(t) is given by series (|3T| ). In the inverse series 

N 

t(s) = £ b k s 2k+l a 2 \ (34) 

k=0 

only terms up to N = 1 are significant to ensure the given accuracy. Coefficients 
b and b\ have the following numerical values 

Consequently, the approximate formula for Fr(s) takes the form 

In the similar way we compute Gr(w), where w = xi(t) — t. The inverse series 
is given by the formula 

N 

. ^ 2k + l 2k-2 

t(w) = 2^ c k w 3 a 3 . (36) 

fc=0 

In this case we have to compute coefficient c k up to iV = 7. We will not present 
here their numerical values. Function Gr(w) takes the form 
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2fc+4 2fc-2 
3 

) 



k=0 



where the coefficients g k have the following approximated numerical values: 

g = 1.1900, g x = 0.0593, g 2 = -0.0299, g 3 = -0.0429, 



g 4 = 0.0037, g 5 = -0.0134, g 6 = 0.0022, g 7 = 0.0002. 

It is important to note that the term k = in Gr(w) has a singular dependence 
on a, i.e. it is proportional to a" 2 / 3 . Such behaviour is caused by the fact that 
term proportional to a , linear in t, is cancelled in the definition of the variable w 
which is proportional to a 2 t 3 for t <C 1. In fact, the singular term cr~ 2 / 3 appears 
already in the series (1361) . 
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3.4 The partial solution f^(x, t) 

For x < t the self-similar solution (1221) is equal to zero, therefore eq. (1291 takes 
the form 

($-^)/f(z,t)-l = 0. (37) 
The general solution of (1371) is of the form 

f I l(x,t)=F L (x + t)+G L (x-t)- l -(x 2 -t 2 ). (38) 

The arbitrary functions F L , Gl and unknown function x (t) can be obtained after 
imposing the following matching conditions: 

f^(t,t) = f*(t,t), f-(x o ,t) = 0, d x f^(x,t)\ x=xo = 0. (39) 

The first of conditions (l39l) gives equality of values of the partial solutions f 1 } and 
/f at the light cone x = t. We do not require equality of spatial derivatives of 
these partial solutions but it turns out that equality of values entails equality of 
derivatives as well. The matching condition at x = t gives an equality F L (2t) + 
Gl(0) = F R (2t), which allows us to obtain the function F L . Without loosing of 
generality we can fix G L (0) = 0, because G L (0) is cancelled in the combination 
F L (x + t) + Gl(x — t). Last two matching conditions (1391 allow us to obtain 
derivatives of functions F L and Gl. Differentiating the second condition in ( |39|) 
with respect to variable t and combining with the third one we obtain 

D+f L {x, t) \ x=s _ t(sl t=t{s) = 0, D^f L (x, t) \ x=w+t{wl t=t{w) = 0, (40) 

where s = x (t) + t, w = x (t) — t and D± = \{d x ±d t ). The solutions of (|40l) 
take the form 

F L (s) = i(a - 2t(s)), G' L (w) = \{w + 2t(w)). (41) 

In order to obtain x (t), we use the equality F L (s) = F' R (s), what gives the 
equation 

a 4 r a 2 , 

— s 5 s 3 - 8s + 16t = 0. 

45 2 

A solution of this equation s(t) can be obtained in the series form. Finally, it gives 
x (t) in the form 

xo(t) = t - itV + ^tV + 0(a e ). (42) 
16 




Figure 4: Trajectory of Xo(t). 



A leading term of deceleration of zero xq is proportional to a 2 for t <^ 1. Fig. 
@] depicts the trajectory of xo(t). The analytical curve which is given by the first 
three terms in formula (1421) is a good approximation for a numerical trajectory if 
t is not greater than t « 2. It gives the limitation for validity of the approximated 
partial solution t). 

The explicit formula for x (t) enables us to obtain the function Gl and the 
solution f^(x, t). Function xo(t), given by (|42|) . is known up to 0(a 4 ) (including 
this term) so t(w) can be obtained up to O(a ). The second formula in (I4TT) gives 



G L (w) 



3 /2\ 1 / 3 



it) 



4/3 _ 



-W 



8 VtrV 135 
Finally, we obtain the partial solution f^(x, t) of the form 



0(a 2 / 3 ). 



(43) 



/L(x ' t) = "^ (151x-119t)(x_t) 



256 



a 



'x + t) 2 + (x + t) 

1 8640 V ; 



3/2 



8 \a 



1/3 



[X 



(44) 
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We note that the singular term ~ <x~ 2 / 3 is present in the formulae (1431) and 
(|44|) . There is nothing unexpected in this fact because for small times the func- 
tion x (t) — t has similar behaviour to the function £i(t) — t from the previous 
paragraph. Some snapshots of the solution f(x, t) for first stage of evolution are 
presented in Figs. [5] and [6] There is a very good agreement between numerical 
and analytical solutions untill t w 2. 
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3.5 Behaviour for later times t > t* 



In the current section we study a numerical solution for the s-K-G model with 
A = —0.1. The solution at t = obeys following conditions 



f(x,0) = -x 2 Q(x), 



d t f(x,t)\ t=0 = 0. 



We focus on times t > t*, where t* « 4.9673. The trajectories of zeros of f(x, t) 
up to t = 25 are depicted in Fig. U\ 
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Figure 7: Trajectories of zeros. 



Zero x (t) decelerates up to time t ~ t*, then it accelerates up to t ~ 2t* when 
it reaches the velocity v — 1. At this moment a very interesting phenomenon 
occurs. Zero x (t) splits into a pair of zeros xo(t) and X\(t) that move similarly to 
x (t) and Xi(t). Numerical values of Xi(t) increase rapidly for t — > 3t*— , which 
suggests that function x\{t) has a vertical asymptote at t = 3t*. This phenomenon 
seems to have an almost periodic character. 

One can show by an elementary calculation (we skip the proof) that the point 
Xo, at which the trivial solution (/ = 0) and a nonzero solution are matched so 
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that the first right-hand spatial derivative of / at x is equal to zero, can not move 
with velocity v > 1. The solution f(x,t), that consists of the partial solutions 
matched at the points that move with the velocity ti / 1, has to be smooth (i.e 
first derivative is continuous) at the matching points. We can conclude that an 
accelerating zero, that is a matching point of the trivial partial solution / = and 
other nonzero partial solution, cannot move faster than v — 1 without changing 
its character. This change means that the first spatial derivative of the field f(x, t) 
is nonzero at this point. It is possible provided that an additional zero, that moves 
with velocity v < 1, appears. The segment of x axis between those two zeros 
widens. Such segment is a support of a partial solution of the second kind. At the 
moment when x reaches the velocity v = 1 the solution f(x, t) returns with good 
approximation to its self-similar initial data (|45l ) (see Fig. [8]). 




Figure 8: Absolute value of the solution f(x, t) at t = 2t* - the solid line, and the 
parabola of self-similar initial data - dashed line. A discrepancy at the right-hand 
side is a pure numerical effect caused by finite size of the grid. 



4 Summary 

A study of perturbed field-theoretic models with V-shaped potential and the scal- 
ing symmetry can give some knowledge that are properties of more general field 
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models with potentials that are sharp at its minima. One of the simplest example 
of such models is the s-G model with perturbation that introduces an additional 
linear term to the field equation, which gives the s-K-G model. We have found 
that such a perturbation is responsible for some new effects. For the same initial 
data the differences between solutions in the s-G and the s-K-G models are caused 
by this linear term. 

The first observation comes from analysis of trajectories of zeros. The zeros of 
the self-similar solutions in the s-G model move with constant velocities, whereas 
velocity of zeros of solutions in the s-K-G model is not monotonous. Moreover, 
their accelerations and decelerations also depend on time. Nevertheless, for small 
times t <C 1 both the acceleration and the deceleration of zeros are proportional 
to the perturbation parameter A. This fact can be useful as a phenomenological 
criterion that allows us to calculate the parameter A from experimental data. 

We can also point out the second qualitative difference between these models. 
There are partial solutions in the s-K-G model that supports expand from zero to 
infinite size within a finite time. It is possible provided that the other partial solu- 
tions disappear simultaneously . For instance, the partial solution <p 1 = j(x 2 — t 2 ) 
at [t, oo) in the s-G model disappears when t — > oo, whereas the partial solution 
/ + in the s-K-G model disappears when t — > ^. 

One of the most interesting results that have been obtained from approximated 
formulae for partial solutions f R (x,t) and ft{x,t) is an observation that these 
partial solutions contain singular terms proportional to A -1 / 3 . It suggests that 
(apart from technical obstacles) the partial solution /_ (x, t) can not be represented 
in the series form in the similar way to f+(x, t). 

The last observation is mainly numerical. The trajectories of zeros of the 
solution f(x,t) are almost periodic. We have proposed here a hypothesis that 
the period has the value 2t* = -, where t* has been obtained from analytical 
calculations - it is the characteristic time for which the trajectory of zero x\ (t) has 
a vertical asymptote. Our hypothesis agrees quite well with the numerical data. 
The periodicity for longer than investigated times is an open question. In our study 
the solutions f(x,0) and f(x, 2t*) are very similar which means that the solution 
f(x, t) returns to the self-similar initial data even though the scaling symmetry is 
broken. 

An important and open question is the behaviour of solutions of the s-K-G 
model for other self- similar initial data or more general initial data. In the group 
of general initial data the most interesting are these for which an initial field con- 
figuration has a finite energy (the energy of the self- similar solutions has an infi- 
nite value). Finally, there are, of course, a variety of perturbations of the potential 
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V(4>) = |0| that can be studied, nevertheless, it is clear that for most of them 
analytical results can be obtained only in approximation. 

5 Acknowledgements 

The author is grateful to Tomasz Tyranowski for his assistance in numerical work 
and valuable discussions as well as to Henryk Arodz, Andrzej Wereszczyhski and 
Jakub Lis for discussions and remarks. 

References 

[1] H. Arodz, P. Klimas and T. Tyranowski, Acta Phys. Pol. B36, 3861 (2005). 

[2] H. Arodz, P. Klimas and T. Tyranowski, Phys. Rev. E73, 046609 (2006). 

[3] PKlimas Acta Phys. Pol. B38, 21 (2007). 

[4] H. Arodz Acta Phys. Pol, B33, 1241 (2002). 

[5] P. Rosenau, J. M. Hyman Phys. Rev. Lett., 70, 564 (1993). 

[6] F. Cooper, H. Shepard and P. Sodano Phys. Rev., E48, 4027 (1993). 

[7] P. Rosenau and A. Pikovsky Phys. Rev. Lett., 94, 174102 (2005). 

[8] C. Adam, J. Sanchez-Guillen and A. Wereszczyhski hep-th/0705.3554 

[9] H. Arodz, P. Klimas and T. Tyranowski preprint hep-th/07 10.2244 

[10] H. Arodz, P. Klimas and T. Tyranowski Acta Phys. Pol, B38, 3099 (2007). 

[11] C. Adam, N. Grandi, J. Sanchez- Guillen and A. Wereszczyhski preprint 
hep-th/07 11.3550 



22 



